Wave localization in strongly nonlinear Hertzian chains with mass defect 
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We investigate the dynamical response of a mass defect in a one- dimensional non-loaded horizontal 
chain of identical spheres which interact via the nonlinear Hertz potential. Our experiments show 
that the interaction of a solitary wave with a light intruder excites localized mode. In agreement 
with dimensional analysis, we find that the frequency of localized oscillations exceeds the incident 
wave frequency spectrum and nonlinearly depends on the size of the intruder and on the incident 
wave strength. The absence of tensile stress between grains allows some gaps to open, which in turn 
induce a significant enhancement of the oscillations amplitude. We performed numerical simulations 
that precisely describe our observations without any adjusting parameters. 
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Wave localization has been studied for long in linear 
systems [l[ or lattices Q and became widely studied more 
recently in nonlinear systems Q. For example, the pres- 
ence of an isotope in a perfect linear crystal is known 
to enhance optical waves absorption at given frequen- 
cies 0] ■ One-dimensional chains of beads interacting via 
the Hertz potential are systems suitable to observe non- 
linear localization effects. A loaded chain of identical 
beads is dispersive, allowing small perturbations to prop- 
agate as linear or weakly nonlinear acoustic waves [5j . In 
contrast, when grains in a chain barely touch one an- 
other, the energy of an impulse only propagates as fully 
nonlinear solitary waves p, 0, 0] resulting from the bal- 
ance between dispersion and nonlinearity of the medium. 
Ncstcrenko early described this regime as a sonic vacuum 
limit [8j . Dissipative effects such as viscoelasticity or fric- 
tion only attenuate and spread these solitary waves @, • 
In contrast, any heterogeneity of the medium capable of 
unbalancing dispersion and nonlinearity results in break- 
ing the solitary wave symmetry. For example, a narrow 
pulse propagating in a chain of beads with decreasing 
sizes develops a long tail which spreads in time the mo- 
mentum transfer [lfj, [ill, Il2 |. Designing powerfull im- 



pact protection systems takes advantage of these fea- 
tures fiol [ill IHi ] . Granular chains made of successions 
of heavy and light beads also proved valuable efficiency 
in energy absorption [Tij ]. More recently, fully nonlinear 
waves with finite- width were observed in chains contain- 
ing periodic mass defects or soft inclusions [15J. Such 
nonlinear dimer chains are expected to support addition- 
nal optical modes and forbidden band gap when sub- 
jected to a static load jl5j |. 

The elementary interaction of either lighter or heavier 
intruders with solitary waves in non-loaded monodisperse 
chains of beads has been investigated numerically [l6j |. 
When a solitary wave reaches a mass defect, energy is 
partially reflected into a backward traveling solitary wave 
and is partially transmitted to the intruder. A heavy im- 



purity slowly translates, leading to a large transmitted 
solitary waves train in the forward direction [lp, simi- 
larly to what was observed in stepped chains [a B 
A light intruder oscillates and scatters forward and back- 
ward weak delayed solitary waves trains [l6| . 

In this letter, we investigate experimentally the inter- 
action of a solitary wave with a light intruder in a non- 
loaded monodisperse chain of beads. We observe that 
the impurity oscillates while interacting with a solitary 
wave. Oscillations remain localized and fade as soon as 
the wave travels away. The presence of a small distance 
gap between the intruder and its nearest neighbors en- 
hances the oscillation amplitude. A multiscale analysis 
predicts the oscillations frequency as a function of in- 
truder size and wave strength. High resolution numerical 
simulations capture all the experimental features. 



Beads are 100C6 steel roll bearings [17j with density 
p = 7780 kg/m 3 , Young's modulus Y = 203 ± 4 GPa 
and Poisson ratio v = 0.3. The chain is made of 20 
equal beads (radius R = 13 mm) and contains an in- 
truder (2.5 mm < Ri < 10 mm) in the middle, as shown 
in Fig. [1] The beads, barely touching one another, are 
aligned on a horizontal Plexiglas track. A special short 
track automatically aligns the intruder on the axis of 
the chain ll|. The chain is ended by a flat, fixed and 
heavy piece of steel. A nonlinear compressive wave (see 
below) is initiated from the impact of a small striker 
(R s = 4 mm). The pulse is monitored by measuring the 
load with a piezoelectric transducer (PCB 200B02, sen- 
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FIG. 1: Upper panel: Experimental setup showing a chain of 
beads with an intruder, sensors and acquisition facilities. 
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sitivity 11.24 mV/N and stiffness 1.9 kN///m) inserted 
inside a bead cut in two parts. The total mass of the 
active bead matches the mass of an original bead and 
the rigidity of the sensor is comparable to bead's mate- 
rial properties. The embedded sensor thus allows non- 
intrusive measurements of the force inside the chain Q. 

Typical force signals measured at the intruder contact, 
plotted in Fig. [2^,1 , shows that large and well defined os- 
cillations appear in the tail of the incident solitary wave. 
The tail corresponds to a slight reflection of the incident 
pulse on the intruder mass defect. Measured oscillations 
in the tail of the force are displayed in more details in 
Figs. [2Jl3-|2^l6 for increasing amplitudes of the incident 
solitary wave. Dissipative processes being negligible pj, 
conservative numerical simulations (see below) shown in 
Fig. [2}d1 demonstrate satisfactory agreements with ex- 
periments, without any adjusting parameter. Calculated 
overlaps between the intruder and neighbor beads, shown 
in Fig. [2)33 (see also the relative displacements in Fig.H^), 
demonstrate that force oscillations correspond to local- 
ized oscillations of the intruder. Figs. [2^2 and[2jD2 depict 
power spectral densities of experimental and calculated 
forces, respectively, in which the high frequency compo- 
nent corresponds to the observed oscillations. 

Fig. [3] presents a longer force acquisition, in which the 
incident pulse and the pulse reflected by the rigid end 
are shown. Force signals exhibit features at two different 
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FIG. 2: Insets (al-a6) are experiments and (bl-b3) are nu- 
merical simulations, both performed in a monodisperse chain 
containing a 3 mm in radius intruder. Dashed line in (al,bl) 
is the 6 beads far the intruder incident solitary wave force 
versus time, and solid line is the force versus time at the left 
intruder's contact. Power spectral densities of these forces are 
plotted in (a2,b2), respectively, showing high frequency con- 
tent in the intruder force. Closer views of oscillating tail in the 
intruder's force are shown in (a3-a6) for increasing incident 
force magnitude (14.1, 17.3, 22.7 and 26.3 N, respectively). 
Solid and dashed lines in (b3) represent overlaps between the 
intruder and left and right neighbor beads, respectively. 
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FIG. 3: Force versus time as detected by the embedded sensor 
in contact with a 3 mm in radius intruder. The first large peak 
in (a) is the incident solitary wave, and the second one is the 
wave after being reflected by the rigid end. The gap opening 
between intruder and nearest beads is indicated in (b) by 
jumps from upper to lower level in the contact switch. High 
and low frequency components of the force, shown in (c), are 
obtained by using a tenth order Butterworth filter with zero 
phase distortion. 

time scales, that are separated by filtering low or high 
frequency contents, as shown in Fig. [3b. Oscillations are 
only observed when the intruder is loaded by the solitary 
wave. In Fig. [3^,, intruder oscillations are only visible 
in the tail of the incident pulse (0.0 ms < t < 1.0 ms). 
Stronger and relatively faster oscillations can be observed 
in the reflected pulse (1.5 ms < t < 3.0 ms), both dur- 
ing the main compression and in the tail. Gap openings 
between the intruder and neighbor beads might explain 
the differences between incident and reflected pulses fea- 
tures. We designed an electrical switch (smooth metallic 
brushes in contact with neighbor beads, connected to a 
fast TTL electrical circuit, see Fig. [IJ to detect any loss 
of contact. The switch is on (1) when beads are in con- 
tact and off (0) when loss of contact occurs. As shown in 
Fig. [3b>, gap first opens right after the incident compres- 
sion (t ~ 0.5 ms and t ~ 2.1 ms) and second at the end 
of tail (f ~ 1.1 ms and t ~ 2.7 ms). The first opening 
is due to the weak reflection of the pulse on the mass 
defect and the second one is due to the momentum and 
energy transfer from the oscillating intruder to neigh- 
bors. Higher amplitude oscillations appears provided a 
gap between the intruder and neighbor exists. 

These observations are corroborated by numerical sim- 
ulation. Calculated displacements of the intruder and 
neighbors shown in Fig. |5l and the overall gap around 
the intruder (the difference between left and right neigh- 
bor beads displacements) shown in Fig. 2b> demonstrate 
that gaps open twice at the intruder, consistently with 
our experiments. Fig. [4b, showing the intruder displace- 
ment relative to the center of mass of the two neigh- 
bors beads, reveals that intruder oscillations even ex- 
ist during the main compression of the incident pulse 
(0.5 ms < t < 1.0 ms, see magnified view in Fig. 0J1). 
Momentum transfer is enhanced by larger strain gradi- 
ent in the presence of gaps. 

Next, we analyzes how oscillations depend on wave 
strength and on intruder size. Frequency of oscillations 
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FIG. 4: Numerical simulations of the incident pulse, under 
same conditions as in Fig.fSbl. (a) Positions versus time of the 
intruder and left and right neighbor beads, (b) Overall gap 
around the intruder demonstrating gap openings occur when 
a solitary wave crosses the intruder, (c) Intruder displacement 
relative to the center of mass of the two neighbor beads, (d) 
Magnified relative displacement during the compression. 



is obtained from the analysis of power spectral densities 
of the force. Following measurements are repeated nine 
times and averaged to minimize errors. We first run a set 
of experiments at constant intruder size (Ri = 2.5 mm) 
while varying the amplitude of the incident solitary wave. 
The force amplitude is obtained from the low-pass filtered 
signals and results are presented in Fig.[5H- We then keep 
the incident pulse amplitude constant (8.6 ± 0.5 N) and 
test several intruder sizes, as depicted in Fig. [SJ}. Exper- 
iments show that the frequency of the intruder nonlin- 
early increases with the incident solitary wave strength 
and depends on the intruder size. 

The physical behavior of solitary waves in chains of 
equal beads and implications for the existence of local- 
ized modes is summarized here. Under elastic defor- 
mation, the energy stored at the contact between two 
elastic bodies submitted to an axial compression cor- 
responds to Hertz potential, Uh — (2/5)kS 5 ^ 2 , where 
6 is the overlap deformation between bodies, k^ 1 — 
(9 + 6')(R- 1 + R'- 1 ) 1 ' 2 and 6 = 3(1 - v 2 )/(4Y) are 
constants, and where R and R' are respective radii of 
curvature at the contact. The force felt at the interface 

3 /2 

derives from Hertz potential, Fjj = dgUn = . In- 
dex + indicates that Hertz force is zero when the beads 
loose contact (no tensile force) . The dynamics of a chain 
of beads is thus described by the following system of N 
coupled nonlinear equations, 



mu ri 



x3/2 
U n ) + 



(tin 



\3/2 
-lJ + 



(1) 



where m and u n are the mass and the position of bead 
n, respectively. Considering long-wavelength perturba- 
tions, such that the strain ip = (—d x u) ~ (6/2R) ~ 
(m/A) -C 1, a continuous equation can be derived from 
Eq. [TJ which admits an exact solution [8j in the form of a 
purely compressive and periodic traveling wave, Tp(x, t) — 
ipm cos 4 [(a; — vt) / (Ry/lO)]. Wave speed, v oc ipU 4 , nonlin- 
early depends on maximum strain. Infinitely small per- 



turbations in the acoustic limit propagate at zero speed, 
acoustic modes are thus forbidden. One hump of this pe- 
riodic function represents a solitary wave solution 0] . In 
addition to analytical estimations, we solve the nonlin- 
ear system of Eqs.[T]by using a fourth order Runge-Kutta 
numerical scheme [18j, the embedded force sensor being 
incorporated in simulations for closer comparisons [l9j |. 
Numerical time step is few orders of magnitude smaller 
than the shortest physical duration and energy conserva- 
tion is fulfilled within a relative error better than 10 . 

The characteristic frequency of localized oscillations 
can be estimated through a multiscale analysis of Eq. [1] 
Here, index n denotes the intruder with radius Ri and 
mass nii — (4/3)tv pR%, and Ki = k(R,Ri) = n(Ri,R) is 
the elastic constant that depends on radii and properties 
of the beads. Index nil denote the two neighbor beads 
with radius R > Ri. We consider two distinct time scales: 
a slow variation of the order of solitary wave duration, 
and a fast time scale associated with intruder oscilla- 
tions. Displacements of the intruder and neighbor beads 
are written as u n ~ u n + u n e lu)t and u n ±i ~ u n ±i, where 
u n and u n are slowly varying functions of time. Harmonic 
oscillations amplitude is assumed negligible compared to 
solitary wave amplitude, u n -c u n . Replacements into 
Eq. Q] leads to, 
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where the first equation, at leading order, provides the 
displacement of the intruder at slow time scale. The sec- 
ond equation, at next order, determines the angular fre- 
quency of the oscillating intruder. We then introduce 
slowly varying forces at the contacts of the intruder, 

3/2 3/2 

F- = Ki(u n -i — u n )+ and F + = Ki(u n -u n +i)+ ■ 
Noticing that they almost behave in phase (F + ~ F_, 
see Fig.[2}33), the first equation indicates that intruder ac- 
celeration oscillates around equilibrium position, u n ~ 0. 
The second equation provides the maximum oscillation 
frequency f m = max[ctj/27r] as a function of the am- 
plitude of the slow time force at the intruder contact, 
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where C = (3/4 7 r v ^rp)/(26»i? 4 ) 1 / 3 . Eq. H shows that 
oscillation frequency tends to zero when the load van- 
ishes: oscillations stops as soon as the solitary wave leaves 
the intruder. Using the characteristics of our beads, 
we find a theoretical estimation, Ct = 2640 Hz/N 1 / 6 . 
Matching Eq. |4] to experiments and simulations shown 
in Fig. we find C e = 2510 ± 151 Hz/N 1 / 6 and C n = 
2531 ±52 Hz/N 1 / 6 , respectively. The agreement between 
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FIG. 5: Inset (a) shows frequency versus force's envelope am- 
plitude, for an Ri — 2.5 mm in radius intruder. Triangles 
and circles correspond to experimental frequency detected in 
the tail of the incident pulse and in the reflected pulse, re- 
spectively. Inset (b) shows frequency versus intruder's radius 
measured in the reflected signal, whose envelope amplitude 
was fixed to F m = 8.7 ± 0.6 N. Dashed lines correspond to 
numerical simulations and straight lines corresponds to Eq.f4] 

experiments, simulations and theory shown in Fig. [5] is 
satisfactory, considering that no adjustable parameters is 
used between experiments and simulations. 

It should be pointed out that weak delayed solitary 
wave trains, responsible for nonlinear leak of oscillating 
energy from the intruder to surrounding (l6j , were de- 
tected in long duration acquisitions of the force recorded 
few beads before or after the intruder (not shown in this 
letter). This attenuation mechanism should lead to an 
exponential decay of the oscillating amplitude within a 
characteristic time much longer than the duration of the 
low frequency force envelope, Tieak/T^p ~ \/ m / TOi S> 1, 
provided the mass of the intruder is lighter than a bead 
of the chain, rm <C m. 

Eq. 2] can also be obtained by considering the low fre- 
quency force envelope amplitude, of the order of F m , 
as a static load at the time scale of fast oscillations. 
The dynamics of a loaded monodisperse chain of beads 
of mass m behaves according to linearized Hertz poten- 
tial around static equilibrium; the intergrain stiffness 

1/3 

reads fc ~ (3/2)k 2 / 3 F to . Small perturbations propa- 
gate as acoustic waves according to the dispersion re- 
lation uj = (2yfc/m)| sin (qR)\. Forcing a single parti- 
cle to move at an angular frequency u> above the cut- 
off uj c = 2y // k/m generates an evanescent wave since 
the wave number is qR = tt/2 — j coshT 1 (u)/u> c ). The 
group velocity v g = (d$t[q]/ 'duj)^ 1 is zero and the per- 
turbation remains localized within a characteristic dis- 
tance x c = = Rj cosh -1 (uj/uj c ). Localization 
achieves when replacing a single particle by a light in- 
truder with mass m, <C m. Roughly estimating the an- 
gular frequency of the intruder from the free oscillation 
frequency between two heavy non-oscillating neighbors, 

LUi ~ y / 2k i /rm where fcj ~ (3/2)/t^ / ' 3 F^ i 3 , leads to Eq. 0] 
expression. The frequency of oscillations exceeds the cut- 
off, oji 3> uj c and the energy remains localized within a 



characteristic distance that does not depend on static 
load and which is smaller than a single radius, x c <C R. 

In conclusion, we have reported observations of energy 
localization in strongly nonlinear chains of elastic beads. 
Localization occurs when a light intruder oscillates at a 
time scale much shorter than the duration of the incident 
solitary wave. Oscillations are excited by the local strain 
gradient and enhanced by the presence of gap at the in- 
truder contact. Oscillations are unable to radiate linear 
acoustic waves as they are forbidden in the sonic vac- 
uum limit. Such nonlinear localization effects are likely 
to appear in polydisperse 3D granular assemblies. For 
instance, a significant number of weakly loaded light par- 
ticles might play an important role on sound trapping at 
high frequency. 
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